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Abstract 

We verify the Kosterlitz Thouless scenario for three different SOS 
(solid-on-solid) models, including the dual transforms of XY-models with 
Villain and with cosine action. The method is based on a matching of 
the renormalization group (RG) flow of the candidate models with the 
flow of a bona fide KT model, the exactly solvable BCSOS model. We 
obtain high precision estimates for the critical couplings and other non- 
universal quantities. 



*Talk presented by K. Pinn at the International Symposium for Lattice Field Theory Tsukuba, 
Japan, November 1991 



For a large class of two-dimensional statistical models, the unambigious verification 
of the Kosterlitz Thouless (KT) scenario is still an open problem. The most recent 
large scale Monte Carlo studies of the XY model |2| clearly favor a KT against a 
second order transition. However, systematic errors are not yet under control. Here 
we present an alternative approach. Our method exploits the fact that the BCSOS 
(body centered solid-on-solid) model can be solved exactly and has been proven to 
exhibit a KT transition ||. The idea is to verify the KT scenario for an SOS model 
by demonstrating that its asymptotic RG flow at criticality matches with the flow 
of the critical BCSOS model. 

The models. References for reviews on SOS models can be found in j|. The 
models to be defined below live on two-dimensional square lattices with periodic 
boundary conditions. The discrete Gaussian model (DG) is dual to the XY model 
with Villain action. The spins hi take integer values. The partition function is 

Z = ^M-K DG £ (fc -/*)*) , 

h <i,j> 

where % and j are nearest neighbor points. The dual of the XY model with cosine 
action also has integer valued spins hi, and partition function 

z=t, n ^i(^ y ), 

h <i,j> 

where the /„ are modified Bessel functions. The ASOS model is defined by the 
partion function 

Z = J2eM-K ASOS £ \hi~hjl). 

h <i,j> 

The BCSOS model was introduced by van Beijeren ||. The lattice is divided in 
even and odd sites, like a checker board. Spins on odd sites take values of the form 
2n + 1/2, spins on even sites are of the form 2n — 1/2, n integer. The partition 
function can be written as 

Z = Y,eM-K BCSOS Y.\h*-h k \), 

h [i,k] 

where i and k are next to nearest (i.e. diagonal) neighbors. Nearest neighbor spins 
hi and hj obey the constraint \hi — hj\ = 1. The critical (roughening) coupling is 
Kr CSOS = \ In 2. The BCSOS model is equivalent to a special case of the six vertex 
model. The critical behavior of several quantities is exactly known and follows the 
predictions of KT theory ||. 

Matching. Universality, first introduced as the coincidence of the critical indices, 
can be understood as a convergence of the RG flow to a universal (model inde- 
pendent) flow as K — > K c and the number of block spin transformations goes to 
infinity. The RG flow can be monitored by studying correlations functions of block 
spins (pi = B~ 2 J2jei hj with increasing block size B. 
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We simulated the BCSOS model at the roughening coupling K^ csos on L x L 
lattices, with L = 12, 16, 24, 32, 48, 64, 96, 128. We measured block spin functions on 
block systems of size I x I, with I — 1,2,4 (L — IB). The statistics was typically a 
few million single cluster updates (see the brief algorithm discussion at the end). 

Motivated by KT theory, we chose as a monitor for the flow of the kinetic term 
the quantities A n = ((<fii — <pj) 2 ), n = 1,2, where i and j are nearest neighbors 
on the block lattice for n — 1, and next to nearest neighbors for n = 2. As a 
monitor for the "fugacity" (periodic perturbation of a massless free field) we chose 
A n+2 = (cos(27m0j)), n = 1, 2, 3 . 

There are two parameters to be adjusted in order to match the RG flow of one 
of the SOS models with that of the critical BCSOS model. First, matching can only 
occur if K sos = K s R os (the roughening coupling of the SOS model). Secondly, one 
has to adjust the ratio of the block sizes b^ s = B sos / B BCSOS . The freedom to set 
^m° 5 7^ 1 is necessary in order to compensate for the different positions of the "bare" 
actions in the Kosterlitz-Thouless flow diagram. Matching occurs if the following 
condition holds: There exists a b^ s and a K^ csos such that for all i and for all I 

aSOS/uSOS -nBCSOS t^SOS\ _^ aBCSOS i tdBCSOS t^BCSOS\ 
A i,l [ b m B ) K R ) ~* A i,l K B i K R ) 

in the limit of large B BCSOS (moderate B's are sufficient in practice). Here Aff s (B, K) 
denotes the observable Aj evaluated in the SOS model at coupling K on a block sys- 
tem of size / x /, each block being of size B x B. Notice that no additional wave 
function renormalization factors are expected, since the minima of the effective po- 
tential of the blocked system are fixed. This is due to the fact that, at all blocking 
levels, the models maintain the global symmetry of shifting all spins by an integer. 

Finite size effects are exactly cancelled since the block systems to be compared 
always have the same number of blocks. Furthermore, since the blocks themselves 
are already large, we expect the matching values of Kff s and b^ s to stabilize for 
small / already. 

In order to determine Kff s we considered, for fixed L sos and /, the following 
two equations: 

Af° s (B sos , Kf os ) = A B f sos (B BCSOS , K BCSOS ) , % = 1, 3 . 

For each of the available values of B BCSOS we solved these equations numerically for 
Kf . The Aff s could be computed for a whole range of couplings with the help 
of the Swendsen-Ferrenberg method ||. Thus we got two values Kf os and K^ os , 
which were in general not identical (matching occurs only for a specific choice of 
b^ s ). In a second step we plotted the values of Kf os and K^ os as function 
of B BCSOS . To obtain continuous curves, we interpolated linearly in \ogB BCSOS . 
The intersection of the two curves K^ os (log B BCSOS ) and K^ os (\ogB BCSOS ) then 
uniquely determined an estimate for the roughening coupling K^ 08 , and, in addi- 
tion, for b^ s . This completes the matching of the SOS and BCSOS flows at the 
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Table 1: (3r and b m = B XY / ' B BCSO for the dual of the XY model with cosine action 
as obtained from the matching of Ai and A 3 



L 


/ 






b 


16 


2 


1.1220(12) 


19.0(1.1) 


0.84(5) 


16 


4 


1.1257(8) 


21.3(4) 


0.75(1) 


24 


2 


1.1214(13) 


26.5(2.3) 


0.91(7) 


24 


4 


1.1225(8) 


28.6(8) 


0.84(2) 


32 


2 


1.1211(12) 


34.5(4.3) 


0.93(10) 


32 


4 


1.1214(8) 


37.4(1.5) 


0.85(3) 


48 


2 


1.1199(11) 


54.1(9.6) 


0.89(13) 


48 


4 


1.1205(7) 


53.6(2.1) 


0.89(3) 


64 


2 


1.1212(11) 


78.(14.) 


0.82(12) 


64 


4 


1.1201(7) 


72.1(3.9) 


0.89(5) 


96 


2 


1.1189(11) 


108.(17.) 


0.89(12) 


96 


4 


1.1194(7) 


100.9(8.0) 


0.95(7) 



roughening transition for given L sos and I. As an example we show in table [I] the 
results for the dual of the XY model with cosine action. 

For all the three SOS models considered, the results for the roughening coupling 
Kr obtained for the various lattice sizes L and sizes / of the blocked system are 
consistent with one another within statistical errors. Only the couplings for I = 4 
on the smallest two lattices sizes and for / = 2 on the smallest lattices slightly 
deviate from the rest. The same is true for the b m . This indicates an extremely fast 
convergence to a universal RG flow of the models. To estimate Kr for the three 
models we averaged the values obtained from the largest L sos we considered both 
for I = 2 and I = 4, and from the second largest L sos for I = 2 only. We arrive at 
the following results: 

(3% Y = 1.1197(5), b* Y = 0.89(5) 
K£ G = 0.6645(6), b° G = 0.31(2) 
Kf os = 0.8061(3), bi sos = 2.8(3) 

The errors are statistical errors. Systematic errors due to deviations from the uni- 
versal flow should be much smaller. 

In order to check the universality of the matching we evaluated the observables 
A 2 , A 4 and A 5 at the critical couplings Kr determined above from Ai and A 3 alone. 
The results show that (within errors) matching occurs for all observables. 

Other non-universal constants. The method allows to determine other non- 
universal constants appearing in the formulas for the divergence of observables near 
the roughening transition 0. For the constants in the asymptotic formula £ = 



4 




K-K R \l/2 
K R I 



) 



we find 



A XY 
A DG 



0.223(13), 

0.078(5), 

0.70(8), 




1.78(2) 
2.44(3) 
1.14(2) 



(for the XY model we used | " ^ | in place of | K K ^ R | )■ 

Our results compare quite well with results from other Monte Carlo simulations 
[0, |J if the systematic errors in these studies are taken into account. Actually, our 
statistical errors are considerably smaller. 

Algorithms. For the algorithms used in the simulation of the DG, the dual XY 
and the ASOS models we used algorithms as described in [|J. Details of the BCSOS 
algorithm, which is also based on the ideas developed in [Q, will be described else- 
where ||. Our BCSOS algorithm has a critical dynamical exponent « 1. Therefore 
more than two thirds of the CPU time we used for this study were spent for the 
BCSOS simulations. With a recently developed new algorithm for the six vertex 
model, the necessary resources would have been considerably smaller ||. 

The matching method was also applied to the Ising interface, for which also a 
cluster algorithm was developed ||10|| . Results on this work will be reported elsewhere 



Acknowledgements. This work was supported in part by the Deutsche For- 
schungsgemeinschaft (DFG), by the German-Israeli Foundation for Research and 
Development (GIF), and by the Basic Research Foundation of The Israel Academy 
of Sciences and Humanities. The numerical simulations were performed at the HLRZ 
in Julich and at the Regionales Hochschulrechenzentrum Kaiserslautern (RHRK). 



5 



References 

[1] W. Janke, K. Nather, Phys. Lett. A157 (1991) 11. 

[2] R. Gupta, C. F. Baillie, preprint LA-UR-91-2364, to appear in Phys. Rev. B. 

[3] R. J. Baxter, Exactly Solved Models in Statistical Mechanics, Academic Press 
(1982). 

[4] H. G. Evertz, M. Hasenbusch, M. Marcu, K. Pinn and S. Solomon, Phys. Lett. 
B254 (1990) 185. 

[5] H. van Beijeren, Phys. Rev. Lett. 38 (1977) 993. 

[6] R. H. Swendsen and A. M. Ferrenberg, Phys. Rev. Lett. 61 (1988) 2635. 
[7] M. Hasenbusch, Ph.D. thesis, Kaiserslautern 1991; 

M. Hasenbusch, M. Marcu, K. Pinn, in preparation. 
[8] M. Hasenbusch, G. Lana, M. Marcu, K. Pinn, in preparation. 

[9] H. G. Evertz, talk presented at this conference; 

H. G. Evertz, G. Lana and M. Marcu, in preparation. 

[10] M. Hasenbusch, S. Meyer, Phys. Rev. Lett. 66 (1991) 530. 

[11] M. Hasenbusch, S. Meyer, in preparation. 



6 



